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Black holes (BHs) in asymptotically anti-de Sitter (AdS) spacetimes have been the subject of 
intense scrutiny, including detailed frequency-domain analysis and full nonlinear evolutions. Re¬ 
markably, studies of linearized perturbations in the time-domain are scarce or non-existing. We 
close this gap by evolving linearized scalar wavepackets in the background of rotating BHs in AdS 
spacetimes. Our results show a number of interesting features. Small BHs in AdS behave as asymp¬ 
totically flat BHs for early/intermediate times, displaying the same ringdown modes and power-law 
tails. As the field bounces back and forth between the horizon and the timelike boundary it “ther- 
malizes” and the modes of AdS settle in. Finally, we have indications that wavepackets in the 
vicinity of fastly spinning BHs grow exponentially in time, signalling a superradiant instability of 
the geometry previously reported through a frequency-domain analysis. 

PACS numbers: 04.70.-s 04.50.Gh 11.25.Wx,11.25.Tq 


I. INTRODUCTION 

The motivation to study BHs in asymptotically AdS 
spacetimes is driven mostly by the gauge/gravity duality. 
This correspondence mia provides a powerful framework 
within which certain strongly coupled gauge theories can 
be mapped to weakly coupled quantum gravitational sys¬ 
tems in one dimension higher. Within this holographic 
setup, a BH is dual to a thermal state. Perturbed BHs 
are dual to non-thermal boundary gauge theories and the 
approach to equilibrium in the gravitational side is trans¬ 
lated to understanding thermalization in the boundary 
gauge theory side mu (see also the review [T2]). 

Understanding the behavior of perturbed BHs in 
asymptotically AdS spacetimes is therefore a central as¬ 
pect of current fundamental and practical research en¬ 
deavors. In asymptotically flat spacetimes, the response 
of a BH spacetime to an external perturbation can be 
roughly divided in three stages. The initial stage corre¬ 
sponds to on-light-cone propagation and is entirely re¬ 
lated to the initial data content; at intermediate times, 
after the initial pulse reaches the null circular geodesic, 
ringdown ensues. In this stage the perturbation decays as 
a series of exponentially damped sinusoids, or quasinor¬ 
mal modes (QNMs) whose properties are intrinsic to the 
BH and do not depend on the external perturbation HU; 
Finally, at late-times backscattering off spacetime curva¬ 
ture becomes important, and a power-law tail becomes 
visible mm- 

It turns out that asymptotically AdS spacetimes in¬ 
troduce novel aspects in the physics of BHs, mostly due 
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to the absence of dissipation at infinity: the boundary 
is timelike and perturbations “bounce off”. Thus, one 
expects the structure of the perturbation spectrum to 
change completely and to even find linear and nonlin¬ 
ear instabilities in these systems [IMS]. A further in¬ 
teresting aspect of the different boundary conditions is 
that, for small AdS BHs (i.e., for BHs whose areal ra¬ 
dius at the horizon is much smaller than the AdS scale) 
one expects imprints of asymptotically flat spacetimes 
in the early signal mm- This result, intuitive from 
a time evolution picture, seems to contradict frequency- 
domain analysis: the mode structure of small AdS BHs is 
completely different from that of say, asymptotically flat 
Schwarzschild or Kerr BHs. Furthermore, understanding 
the linear regime of perturbed BHs is important as bench¬ 
marks for nonlinear evolutions, some of which have been 
recently reported [221424] . Surprisingly, time-domain evo¬ 
lutions in the background of rotating BHs in AdS space- 
times are missing. This work intends to close this gap. 


II. FRAMEWORK AND NUMERICAL 
PROCEDURE 


A. Wave Equation 


The massless Klein-Gordon equation in Kerr-AdS 
background describes scalar field perturbations in the 
space-time of Kerr-AdS BHs mm- In Boyer-Lindquist 



2 


coordinates, this equation takes the form[]] 
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where M is the BH mass, a its angular momentum per 
unit mass, A r = (r 2 + a 2 )(l + js) — 2Mr , Ag = 1 — 

2 n 2 

j 2 cos 6 and £ = 1 — ^ ■ We say that a BH is small when 
M/i < 1. This definition is equivalent to requiring that 
the areal radius at the horizon is small, when compared 
to the AdS radius m- 

The equation above is a linear, hyperbolic, homo¬ 
geneous (3+l)D partial differential equation that was 
shown to yield a well-posed initial-value problem in 
Refs. 2SU2B1- 


B. Evolution Code 


A stable evolution scheme for the Teukolsky master 
equation of Kerr spacetime was developed in Ref. m 
using the well-established, 2-step Lax-Wendroff numeri¬ 
cal evolution scheme. Our evolution code uses the exact 
same approach, therefore the contents of this section are 
largely a review of the work presented in the literature. 

Our code uses the tortoise coordinate r* in the radial 
direction, which is related to the Boyer-Lindquist coor¬ 
dinates r through the transformation 

dr* = ^ dr ■ ( 2 ) 

These coordinates are better suited for performing nu¬ 
merical evolutions in a Kerr space-time background for a 
number of reasons that are detailed in Ref. E3- Next, we 
factor out the azimuthal dependence and use the ansatz, 

V(t,r*,0,4>) = e im<l, ${t,r*,9) (3) 


that allows us to reduce the dimensionality of the orig¬ 
inal (3+l)D Eqn. [l] to a system of (2+l)D equations. 
Defining 
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allows the Teukolsky equation to be rewritten in first 
order form as 

d t u -(- Md r *u + Lu -f Au = 0, (7) 

where 

m = {d)^, $7,n^,n/} (8) 


is the solution vector. The subscripts R and I refer to 
the real and imaginary parts respectively (recall that the 
Teukolsky function U/ is a complex valued quantity). Ex¬ 
plicit forms for the matrices M , A and L can be easily 
computed given the expressions we enlist above. Rewrit¬ 
ing Eq. 0 as 
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S = -(M - D)d r *u- Lu- Au, (11) 


and using the Lax-Wendroff iterative scheme, we obtain 
stable evolutions. Each iteration consists of two steps: 
In the first step, the solution vector between grid points 
is obtained from 

<ti/2= 2 «+!+<)- ( 12 ) 
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This is used to compute the solution vector at the next 
time step, 
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The angular subscripts are dropped in the above equation 
for clarity. All angular derivatives are computed using 
second-order, centered finite difference expressions. 

We set $ and n to zero on the outer radial boundary, 
and “ingoing” on the inner boundary (which is placed 
extremely close to the hole’s horizon i.e. at a large neg¬ 
ative value of the radial r* coordinate). Symmetries of 
the spherical harmonics are used to determine the angu¬ 
lar boundary conditions: For even \m\ modes, we have 
dg ( E> = 0 at 6 = 0,7r while $ = 0 at 8 = 0,7r for modes of 
odd \m\. 


C. Initial data 


1 Notice the minor typo in Eq.(2.6) of Ref. m, the sin 2 0 in the 
denominator of the second term in the r.h.s should be a sin0. 


We specify initial data for the evolution code by simply 
setting to zero everywhere and choosing an ordinary 
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Gaussian radial-profile for II and a pure spherical har¬ 
monic in the angular dimension. More specifically, we 
choose the initial data to be of the form: 

ne'^ = e^Y lm (14) 

so that the Gaussian is centered at r* = 0 and has a 
width A = 4. We always choose l = m in the current 
work. 


D. Parameters 

For the numerical evolutions, we use a radial grid 
that is uniform in the tortoise coordinate r*. The in¬ 
ner boundary is placed at a large negative value of r* 
in a manner that it is outside of the domain for causal 
influence for the location and evolution duration of inter¬ 
est. The outer boundary is located precisely at r* = f £, 
or r = oo in Boyer-Lindquist coordinates. Everything is 
computed in units of BH mass M which is numerically 
chosen to be 1.0. Our numerical computations cover a 
range of parameters: a/M in the range of 0.0 — 0.995; 
t in the range of 20 — 10 4 ; and m with values 0, 1 and 
2. For simplicity we will focus on results i = 64,10 4 
and a/M = 0, 0.995. The numerical code is fully par¬ 
allelized to operate on a Beowulf cluster using domain- 
decomposition along the radial dimension, and exhibits 
excellent scaling. We typically use radial grid resolutions 
of M/40 and angular resolutions of 7r/64, with a Courant 
factor of 0.5 for the time-stepping. Each simulation was 
completed using 500 processor-cores each, taking a few 
days of walltime. 

III. RESULTS 



FIG. 1. Very long time evolution of a spherically symmetric 
(Z = m = 0) wavepacket in a Kerr-AdS background, with 
£ — 10 4 M, contrasted against the evolution of the same initial 
wavepacket in Schwarzschild-AdS. 



FIG. 2. Time evolution of a spherically symmetric (Z = m = 
0) wavepacket in a small Schwarzschild-AdS background, with 
£ = 10 4 iU, contrasted against the evolution of the same initial 
wavepacket in asymptotically flat spacetime. The pulse shows 
signs of Schwarzschild ringdown and power-law tail at very 
early times. The complete agreement with that of asymptoti¬ 
cally BHs is shown in the inset, which zooms-in the evolution 
at early times. At very late times, not shown here, the two 
evolutions differ appreciably. 



FIG. 3. Time evolution of a wavepacket in a small Kerr- 
AdS background, with £ = 64 M, for an initial l = m = 0 
scalar wavepacket. The BH is spinning at a rate a/M = 0 
(black line) and a/M = 0.995 (blue line). The m = 0 pulse 
“thermalizes” very quickly, and at late times, we find good 
agreement with the Kerr-AdS QNMs computed via frequency- 
domain methods m- The inset shows an early-time zoom-in 
of the evolution, which shows signs of ringdown similar to 
that of asymptotically flat BHs; however, before ringdown 
“dies off,” the signal has had time to travel to and back from 
the boundary, significantly distorting the waveform. 


We focus now on fields extracted at the equator 6 = 
7t/ 2 and at r = 20 M. A typical waveform for axially 
symmetric (m = 0) fluctuations is shown in Figs. [l|3] for 
£ = 10 4 M, 64M. The early time evolution is dictated by 
the initial data and by the scattering off the BH hori¬ 
zon. Accordingly, the three stages - prompt response 
followed by ringdown and a power-law tail - of linearized 
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perturbations discussed in the introduction are seen here 
very clearly. In fact, at early times t <C £ and for very 
small M/£, the quantitative details of these stages are 
identical to the ones in asymptotically flat spacetime. In 
particular small BHs in AdS excite the same ringdown 
modes and tails as those in Kerr backgrounds. This is 
clearly seen for small BHs such as the ones in Figs. HE a 
quantitative comparison of the signal with that of asymp¬ 
totically flat spacetimes is provided in Fig. [2j where we 
compare the signal against that of a BH with the same 
mass but living in asymptotically flat spacetime. The two 
signals are identical up to timescales of the order of the 
AdS radius. Such results are not surprising: the QNMs 
of asymptotically flat BHs are modes leaking from the 
unstable null circular geodesic [28] and therefore could 
be expected to be excited as well. Perhaps the unex¬ 
pected feature is that they play no role whatsoever in 
the spectral content of the wave operator in AdS, or in 
other words, the QNMs of asymptotically flat BHs do not 
show up as poles of the retarded Green’s function in AdS 
(see also Ref. :20, ETj for a related phenomena in con¬ 
nection with BHs in astrophysical environments). As is 
well known, a frequency-domain analysis for the QNMs 
of small BHs in AdS returns extremely long-lived modes 
lying close to the purely real, AdS normal modes 

u>£ = 3 + l + 2n, (15) 

with l, n = 0,1, 2,.. angular and radial numbers respec¬ 
tively mm- 

To summarize, imprints of asymptotic flatness are 
present in waveforms of asymptotically AdS BHs. The 
overall magnitude of these imprints can be varied, by ap¬ 
propriately tuning the initial data (for example, it can 
be chosen in such a way that ringdown excitation is sup¬ 
pressed (see Refs. [2§1I5U] for studies in frequency-domain 
and Refs. [3Tl - i34] for corresponding analysis in the time- 
domain), but is otherwise generic. This is one of the main 
messages of our work. 

Because the AdS radius introduces an extra length 
scale in the problem, the evolution is characterized by 
two extra stages relatively to its asymptotically flat coun¬ 
terpart. The first occurs on timescales large enough that 
AdS curvature becomes important, modifying the char¬ 
acteristic of the “tail”. Because the power-law tails are a 
large wavelength effect, they feel the impact of the AdS 
modifications and our results suggest that these changes 
will eventually contribute to the settling in of the AdS 
resonant modes. The final stage is the reflection off the 
boundary at t ~ £ir/2. These two stages are apparent 
already in Fig. |T| where the reflection off the bound¬ 
ary is visible. However, a smooth transition to the AdS 
regime would take a large number of reflections and con¬ 
sequently long timescales. We therefore focus on a some¬ 
what larger AdS BH, with £ = 64M, shown in Fig. [3] for 
a = 0, 0.995 M. Furthermore, the inset in Fig. [2] shows 
that ringdown lasts for ~ 80M; thus, similar initial data 
will result in a ringdown stage similar to that of asymp¬ 
totically flat BHs, as long as l 80 M, i.e., as long as 


the signal does not have time to travel back from the 
boundary. Thus, we expect that with £ = 64 M the ring- 
down stage will not be identical to its asymptotically flat 
counterpart. 

Both these features are borne out of our results, shown 
in Fig. [3j The inset, zooming in on the evolution at 
early times, shows that the “asymptotically-flat” mem¬ 
ory (i.e., ringdown characteristic of asymptotically flat 
BHs) is quickly distorted; the reason, as we explained, 
is that on timescales of ~ 80 M the field bounced back 
from the boundary already. We find that the transition 
to the AdS regime occurs after a few reflections off the 
boundary. Notice that upon each reflection at the bound¬ 
ary the wave is partially absorbed by the BH, leading to 
a steady decrease in amplitude HJ We find that both 
the resonant frequencies and damping times are well ap¬ 
proximated by frequency-domain results [12]; in partic¬ 
ular, for small BHs one recovers (15). For example, a 


Fourier analysis of the m = 0, a = 0 waveform yields 
£uj ~ 2.98 — i0.0045, in very good agreement with nu¬ 
merical and analytical results in the frequency-domain, 
which yield tw ~ 2.92 — j0.0049 pR^, R). There are no 
power-law tails appearing at time scales much larger than 
£ because of the decay properties of the potential SHU- 
The overall picture is the same also for non- 
axisymmetric modes, whose evolution is summarized 
in Figs. 4][5 However, because these modes are non- 
axisymmetric, they are prone to superradiant effects 
when interacting with the BH [ism Eg: for small 
enough rotating BHs in AdS, the resonant frequency is 
that of pure AdS, ui ~ l/£ and can be superradiant, i.e., 


u) < mfl. 


(16) 


with H the horizon’s “angular velocity.” In this regime, 
after reflection at the boundary the wave is amplified close 
to the horizon, leading to an exponential growth of the 
field. This is usually known as a BH bomb or super¬ 
radiant instability PEMEI ESI EH- Figure [5] shows indi¬ 
cation that the evolution of wavepackets around small, 
rapidly rotating BHs in AdS does lead to a growth of 
the field. The initial stages show the same features as 
thhose already discussed, with the initial pulse scatter¬ 
ing off the BH and bouncing of the AdS boundary. This 
process will typically filter high-frequency components off 
the signal, because they are easily absorbed by the BH. 
At intermediate-late stages, the signal consists mainly on 
lower-frequency components which are prone to super¬ 
radiance growth. We estimate the growth-timescale to 
be of order 4 x 10 2 * * * 6 M, which agrees with the estimate of 
Ref. .[35] and with frequency-domain calculations HMg. 


2 In fact, this simple picture of waves at the boundary and con¬ 

sequent absorption at the horizon can be used to estimate the 

decay rate of any fluctuation in a confined geometry solely from 

the knowledge of the absorption probabilities of asymptotically 

flat BHs [35 ] 
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FIG. 4. Time evolution of a dipolar l = 1 wavepacket in a small Kerr-AdS background, with l = 10 4 A7. Left panel shows 
m = 1 waveforms in the background of a BH with spin a = 0 and a = 0.995 M. The right panel zooms in on the corresponding 
signal at late times. The waveform begins with a high frequency part corresponding to a Schwarzschild ringdown. At timescales 
t ~ £ the effects of the cosmological constant are felt and a low frequency component eventually settles in before complete 
reflection at the boundary. Eventually, the low frequency component dominates the signal. The turn-over to the AdS, with 
“thermalization” at low frequencies is more clearly seen in Fig. [3] 



time 

FIG. 5. Very long time evolution of an m = 1 wavepacket 
in a Kerr-AdS background, with t = 64 M and a/M = 0.995, 
contrasted against the evolution of the same, however, with 
m = 0. Instability on the scale of 10 6 M develops in the non- 
axisymmetric case, as predicted by frequency-domain analy¬ 
sis fl5Hl8| . 


IV. CONCLUSIONS 

The early-time evolution of fluctuations around small 
AdS BHs is identical to that of their asymptotically flat 
counterpart: the evolution can be divided in the same 
three stages, prompt response, ringdown and late-time 
power-law tail. We highlight that these stages are not 
captured by a frequency-domain analysis [201 [5T]. What 
frequency-domain analysis does capture are the late-time 


dynamics, when the fluctuations have had time to bounce 
off the boundary several times and eventually decay in 
the Kerr-AdS quasinormal modes. These modes decay in 
time and can be understood as a sequence of reflections 
and subsequent absorption at the BH horizon. Finally, 
we have have indications - directly in the time-domain 
- that small, rotating BHs in AdS are unstable against 
a superradiant mechanism in the sense that small non- 
axisymmetric fluctuations grow nsmansi- 

The dynamics in asymptotic AdS spacetimes are rich 
but complex. Their understanding requires the full non¬ 
linear machinery of Einstein equations, complemented 
with robust linearized analysis. Our study intends to 
partially fill the gap in linearized, time-evolution studies 
in asymptotic AdS spacetimes. 
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